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It has been shown experimentally a long time ago that the magnetic ordering 
causes an anomalous behavior of the electron resistivity in ferromagnetic crystals. 
Phenomenological explanations based on the interaction between itinerant electron 
spins and lattice spins have been suggested to explain these observations. We show 
by extensive Monte Carlo simulation that this behavior is also observed for the 
resistivity of the spin current calculated as a function of temperature (T) from low-T 
ordered phase to high-T paramagnetic phase in a ferromagnet. We show in particular 
that across the critical region, the spin resistivity undergoes a huge peak. The origin 
of this peak is shown to stem from the formation of magnetic domains near the phase 
transition. The behavior of the resistivity obtained here is compared to experiments 
and theories. A good agreement is observed. 

PACS numbers: 

I. INTRODUCTION 

The interplay between magnetic properties and transport behavior has been a subject 
of continuous experimental investigations.-^ It has been shown in these works that the 
resistivity has an anomalous peak in the magnetic transition temperature region. De Gennes 
and Friedel,- and later Fisher and Langer,- have suggested that the peak at the transition is 
due to the magnetic short-range interaction. Various analytical methods such as mean-field 
approximationsS^have been employed to explain experimental observations. 

The problem of spin-dependent transport has been also extensively studied in magnetic 
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thin films and multilayers. The so-called giant magnetoresistance (GMR) was discovered 
experimentally twenty years ago.- 1 ^ Since then, intensive investigations, both experimentally 
and theoretically, have been carried out.— >^ Experimental observations show that when the 
spin of an itinerant electron is parallel to the spins of the environment it will go through 
easily while it will be stopped if it encounters an antiparallel spin medium. The resistance 
is stronger in the latter case resulting in a GMR. Although many theoretical investigations 
have been carried out, to date very few Monte Carlo (MC) simulations have been performed 
regarding the temperature dependence of the dynamics of spins participating in the cur- 
rent. In a recent work,— we have investigated by MC simulations the effects of magnetic 
ordering on the spin current in magnetic multilayers. Our results are in agreement with 
measurements.— 

We study in this paper the transport of itinerant electrons traveling inside a bulk ferro- 
magnetic crystal using extensive MC simulations. The electron spin is supposed to be the 
Ising model. Various interactions are taken into account, in particular interaction between 
itinerant spins and interaction between itinerant spins and lattice spins. 

The paper is organized as follows. Section II is devoted to the description of our model and 
the rules that govern its dynamics. In section III, we describe our MC method and discuss 
the results we obtained. These results are in agreement with experiments and theories. 
Concluding remarks are given in Section IV. 

II. MODEL 

We consider in this paper a bulk ferromagnetic crystal. We use the Ising model and the 
face-centered cubic (FCC) lattice with size 4N X x N y x N z . Periodic boundary conditions 
(PBC) are used in the three directions. Spins localized at FCC lattice sites are called "lattice 
spins" hereafter. They interact with each other through the following Hamiltonian: 

W^-J^-S,-, (1) 

where Si is the Ising spin at lattice site i, Yluj) indicates the sum over every nearest-neighbor 
(NN) spin pair (Sj, Sj), J(> 0) being the NN interaction. 

In order to study the spin transport in the above system, we consider a flow of itiner- 
ant spins interacting with each other and with the lattice spins. The interaction between 
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FIG. 1: Lattice magnetization versus temperature T. T c is ~ 9.75 in unit of J = 1. 
itinerant spins is defined as follows, 

7~£m = ^ t KijSj, ■ Sjj (2) 

where is the Ising spin at position r^, and Y^u,j) denotes a sum over every spin pair (s^ Sj). 
The interaction i£y depends on the distance between the two spins, i.e., r^ = \fi — r}|. 
A specific form of will be chosen below. The interaction between itinerant spins and 
lattice spins is given by 

T~C r = Iid s i ' ^i' (^) 

where the interaction /^j depends on the distance between the itinerant spin Sj and the 
lattice spin Sj. For the sake of simplicity, we assume the same form for Kij and ijj, namely, 
Ki j = K exp(— rij) and I^j = Iq exp(— r^ ), where K and Jo are constants. 

The procedure used in our simulation is described as follows. First we study the thermo- 
dynamic properties of the bulk system alone, i.e., without itinerant spins, using Eq. (CQ). We 
perform MC simulations to determine quantities as the internal energy, the specific heat, 
layer magnetizations, the susceptibility, ... as functions of temperature T.— From these 
physical quantities we determine the critical temperature T c below which the system is in 
the ordered phase. We show in Fig. [T]the lattice magnetization versus T. 

Once the lattice has been equilibrated at T, we inject iVo itinerant spins into the system. 
The itinerant spins move into the system at one end, travel in the x direction, escape 
the system at the other end to reenter again at the first end under PBC Note that PBC 
are used to ensure that the average density of itinerant spins remains constant during the 
time (stationary regime). The dynamics of itinerant spins is governed by the following 



interactions: i) an electric field E is applied in the x direction. Its energy is given by 
He = — E- v, where v is the velocity of the itinerant spin; ii) a chemical potential term which 
depends on the concentration of itinerant spins within a sphere of radius D 2 ("concentration 
gradient" effect). Its form is given by 7i c = Dn(r), where n(r) is the concentration of 
itinerant spins in a sphere of radius D 2 centered at r. D is a constant taken equal to Kq for 
simplicity; iii) interactions between a given itinerant spin and lattice spins inside a sphere of 
radius D\ (Eq. [3]); iv) interactions between a given itinerant spin and other itinerant spins 
inside a sphere of radius D 2 (Eq. 

Let us consider the case without an applied magnetic field. The simulation is carried out 
as follows: at a given T we calculate the energy of an itinerant spin by taking into account 
all the interactions described above. Then we tentatively move the spin under consideration 
to a new position with a step of length vq in an arbitrary direction. Note that this move is 
immediately rejected if the new position is inside a sphere of radius r centered at a lattice 
spin or an itinerant spin. This excluded space emulates the Pauli exclusion principle in the 
one hand, and the interaction with lattice phonons on the other hand. If the new position 
does not lie in a forbidden region of space, then the move is accepted with a probability 
given by the standard Metropolis algorithm.— 

III. MONTE CARLO RESULTS 

We let N Q itinerant spins travel through the system several thousands times until a 
steady state is reached. The parameters we use in most calculations, except otherwise 
stated, s = S = 1 and N x = N y = N z = 20. Other parameters are Dy — D 2 = 1 (in unit 
of the FCC cell length), Kq = Iq = 2, Nq = 20 3 (namely one itinerant spin per FCC unit 
cell), Vq — 1, ro = 0.05. At each T the equilibration time for the lattice spins lies around 10 6 
MC steps per spin and we compute statistical averages over 10 6 MC steps per spin. Taking 
J = 1, we obtain T c ~ 9.75 for the estimate of the critical temperature of the lattice spins 
(see FigJTJ. 

In Fig. [2] we sketch the characteristic traveling length A computed after a fixed lapse of 
time as a function of temperature T. As can be seen, A is very large for T < T c . We note 
that there is a small depression in the transition region. We will show below that this has 
an important consequence on the spin current: the resistance undergoes a cusp. 
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FIG. 2: Characteristic traveling length A in unit of the FCC cell length versus temperature T, for 
100 MC steps. 



FIG. 3: Resistance R in arbitrary unit versus temperature T. Our result using the Boltzmann 
equation is shown by the continuous curve. 

We define the resistance R as R = -, where n is the number of itinerant spins crossing a 
unit area perpendicular to the x direction per unit of time. In Fig. [3] we show the resistance 
R clS 8b function of temperature. 

There are several striking points observed here. First, R is very low in the ordered 
phase and large in the paramagnetic phase. Below the transition temperature, there exists a 
single large cluster with some isolated defects, so that any itinerant spin having the parallel 
orientation goes through the lattice without hindrance. The resistance is thus very small. 
Secondly, R exhibits a cusp at the transition temperature, the existence of which is due to 
the critical fluctuations in the phase transition region. De Gennes and Friede^ and later 
Fisher and Langer- and then Kataoka£ among others have shown that the resistance is 
related to the spin correlation < SiSj > which is nothing but the magnetic susceptibility. 
The latter diverges at the transition. This explains at least partially the cusp of R observed 
here. Another picture to explain the cusp is that when T c is approached large clusters of up 
(resp. down) spins form in the critical region above T c . As a result, the resistance is much 
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larger than in the ordered phase: itinerant electrons have to steer around large clusters in 
order to go through the entire lattice. Thermal fluctuations are not large enough to allow 
the itinerant spin to overcome the energy barrier created by the opposite orientation of the 
clusters in this temperature region. Of course, far above T c , most clusters have a small 
size, so that the resistivity is still quite large with respect to the low-T phase. However, 
R decreases as T is increased because thermal fluctuations are more and more stronger to 
help the itinerant spin overcome energy barriers. Note in addition that the cluster size is 
now comparable with the radius Di of the interaction sphere, which in turns reduces the 
height of potential energy barriers. We have checked this interpretation by first creating an 
artificial structure of alternate clusters of opposite spins and then injecting itinerant spins 
into the structure. We observed that itinerant spins do advance indeed more slowly than 
in the completely disordered phase (high-T paramagnetic phase). We have next calculated 
directly the cluster-size distribution as a function of T (not shown) using the Kopelman 
algorithm.— The result confirms the effect of clusters on the spin conductivity. The reader 
is referred to our previous work— for results of a multilayer case. Note that in Fig. [3] we 
also show the results from our theory using the Boltzmann equation in the relaxation-time 
approximation. The agreement is remarkable. We will show the details of this theory 
elsewhere.— 



IV. CONCLUDING REMARKS 



We have shown in this paper first results of MC simulations on the transport of itinerant 
spins interacting with localized lattice spins in a ferromagnetic FCC crystal. Various inter- 
actions have been taken into account. We found that the spin current is strongly dependent 
on the lattice spin ordering: at low T itinerant spins whose direction is parallel (antiparallel) 
to the lattice spins yield a strong (weak) current. At high temperatures, the lattice spins 
are disordered, the current of itinerant spins is very weak and does not depend on the input 
orientation of itinerant spins. As a consequence, the resistance is very high at high T. We 
found in the transition region between low-T and high-T phases a peak of the resistance 
which is due to the existence of domains of lattice spins. Experiments and theories support 
our finding. 



The authors are grateful to Sylvain Reynal for discussions. 



1 F. C. Shwerer and L. J. Cuddy, Phys. Rev. 2, 1575 (1970). 

2 Alia E. Petrova, E. D. Bauer, Vladimir Krasnorussky, and Sergei M. Stishov, Phys. Rev. B 74, 
092401 (2006). 

3 S. M. Stishov, A.E. Petrova, S. Khasanov, G. Kh. Panova, A.A.Shikov, J. C. Lashley, D. Wu, 
and T. A. Lograsso, Phys. Rev. B 76, 052405 (2007). 

4 P.-G. de Gennes and J. Friedel, J. Phys. Chem Solids 4, 71 (1958). 

5 M. E. Fisher and J.S. Langer, Phys. Rev. Lett. 20, 665 (1968). 

6 Mitsuo Kataoka, Phys. Rev. B 63, 124435 (2001). 134435-1 (2001) 

7 A. L. Wysocki, K. D. Belashchenko, J. P. Velev, and M. van Schilfgaarde, J. Appl. Phys. 101, 
09G506 (2007). 

8 M. N. Baibich, J. M. Broto, A. Fert, F. Nguyen Van Dau, F. Petroff, P. Etienne, G. Creuzet, 
A. Friederich and J. Chazelas, Phys. Rev. Lett. 61, 2472 (1988). 

9 P. Grunberg, R. Schreiber, Y. Pang, M. B. Brodsky and H. Sowers, Phys. Rev. Lett. 57, 2442 
(1986); G. Binash, P. grunberg, F. Saurenbach and W. Zinn, Phys. Rev. B 39, 4828 (1989). 

A. Barthelemy et al, J. Mag. Mag. Mater. 242-245, 68 (2002). 

1 See review by E. Y. Tsymbal and D. G. Pettifor, Solid State Physics (Academic Press, San 
Diego), Vol. 56, pp. 113-237 (2001). 

2 K. Akabli, H. T. Diep and S. Reynal, J. Phys.: Condens. Matter 19, 356204 (2007). 

3 R. Brucas and M. Hanson; J. Magn. Magn. Mater. 310, 2521 (2007). 

4 K. Binder and D. W. Heermann, Monte Carlo Simulation in Statistical Physics, Springer- Verlag, 
New York (1988). 

5 J. Hoshen and R. Kopelman, Phys. Rev. B 14, 3438 (1974). 

6 K. Akabli and H. T. Diep, to be submitted for publication. 



